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INTRODUCTION 


Requirements  for  controlled  structures  in  space  for 
astronomy,  communication  networks,  near-earth  scientific  studies, 
and  space  solar  power  alternatives  have  recently  generated  much 
interest  in  the  active  control  of  flexible  structures  [l].  Such 
applications  typically  call  for  very  large  lattice  structures  (up 
to  several  kilometers  in  extent)  which  would  support  arrays  of 
antennae,  solar  cells,  etc.  A  lattice  is  usually  constructed  by 
assembling  many  simple  truss-like  substructures  into  a  large 
periodic  structure. 

The  scope  of  the  present  study  is  limited  to  the  active 
control  of  the  propagation  of  vibrational  disturbances  in  a 
single  truss-like  substructure,  many  of  which  might  comprise  a 
lattice  structure.  Reduction  of  the  problem  to  this  scale 
greatly  simplifies  the  vibrational  dynamics,  as  the  effects  of 
counaary  reflections  are  omitted.  The  problem  is  further 
simplified  by  idealizing  the  substructures  as  thin  elastic  rods 
wnich  are  subjected  only  to  longitudinal  stress  wave 
disturbances . 

This  study  investigates  the  concept  of  utilizing  a 
feeuforeward  controller  to  sense  and  subsequently  cancel  incoming 
longitudinal  stress  wave  disturbances  by  superposing  the  negative 
of  each  individual  disturbance.  The  objective  cf  such  a 
controller  is  to  isolate  a  portion  of  the  substructure  from 
vibration  due  to  incoming  disturbances.  While  the  report  is 
primarily  aevoted  to  deriving  and  analyzing  the  response 


characteristics  of  a  substructure  with  such  a  feedforeward 


controller  implemented,  a  discussion  is  provided  in  Appendix  A  on 
the  response  of  a  counterpart  feedback  controller.  It  is  r.oped 
that  once  developed,  the  ability  to  control  the  propagation  of 
disturbances  in  single  substructural  elements  may  have  an 
application  to  the  problem  of  lattice-wide  vibration  control. 


and  p  do  not  vary  with  x 


(2)  The  substructure  shape  is  prismatic,  so  that  the  cross 
sectional  area  A  does  not  change  with  x. 

(3)  Uniaxial  stress  is  assumed,  and  although  uniaxial  strain 
is  not  assumed,  the  inertial  effects  associated  with 
lateral  strain  (Poisson  effect)  are  neglected. 


(4;  Stresses  in  the  substructure  are  assumed  to  be  subject  to 
no  attenuation.  This  assumption  is  not  unrealistic, 
since  the  vibration  of  many  long,  flexible,  lightweight 
structures  in  space  is  very  lightly  damped  [l]. 


A  useful  solution  to  eqn .  (3)>  credited  to  D'Alembert 

(1747),  is  [2] 


u(x,t) 


(4) 


where  f  and  f  are  arbitrary  functions  which  represent 

12 

disturbances  propagating  at  a  velocity  c^  ;  f  ^  is  associated  with 
the  components  of  the  total  disturbance  which  are  travelling  in 

the  positive  x  direction,  and  f  represents  the  components 

2 

travelling  in  the  negative  x  direction.  The  velocity  of 

longitudinal  wave  propagation,  c  ,  is  determined  by  the  material 

o 

properties  of  the  elastic  solid,  and  is  given  by  [2] 


Stress  wave  propagation  can  be  shown  to  follow  the  same  rule 
of  propagation.  Taking  the  partial  derivative  of  eqn.  (3)  with 


respect  to  position  gives 


J_  /  32u(x,t)  \  _3_  /  p  32u(x, t)  \ 
3x  \  3x2  /  3x  \  E  3t2  / 


which  may  be  rewritten  as 


32  /  3u(x,  t)  \  _  32  /  p  3u(x,  t)  \ 

3X2  \  3x  /  3t2  \  E  3x  / 


Next,  eqn.  (2)  is  used  to  write 


32e(x,t) 


32£(x,t) 


E  3f 


which,  from  Hooke's  Law  and  the  homogeneity  criterion,  becomes 


32o (x, t)  p  32a(x,t) 


The  D'Alembert  solution  form  thus  pertains  for  stress  wave 
propagation  as  well,  so  that  it  is  possible  to  write 


0<x-t>  ’  8,(t")  +8;(t+t) 


where  g  and  g  are  arbitrary  functions  that  will  be  specifically 

12 

determined  by  the  initial  conditions  or  forcing  function  in  a 
given  problem. 

Lqn .  (10)  describes  stress  waves  which  propagate  without 

distortion  or  attenuation  through  the  substructure.  Thus,  input 
waves  of  arbitrary  "shape"  propagating  in  the  positive  x 
direction  past  point  x  in  the  rod  of  Fig.  1  will  later  pass 

l 

point  x  unchanged  from  their  original  configuration. 

2 

This  nondispersive  characteristic  may  be  expressed  as 


a  (x2,t)  =  a  (x  ^ ,  t-x12) 


where  a  (x2,t)  denotes  the  magnitude  of  a  stress  wave  at  cross 
section  x^  at  time  t,  as  the  wave  propagates  through  the 
substructure.  The  superscript  "  +  "  indicates  that  propagation  is 
in  the  positive  x  direction,  and  t  is  the  time  constant 
associated  with  the  "transportation  lag"  f 3 J  of  propagation  along 

a  length  &  of  tne  substructure  between  sections  x  and  x 

12  12 


x  -  x 
2  1 


Nondispersive ,  unattenuated  propagation  is  illustrated  in 

Fig.  2,  which  shows  a  schematic  of  a  stress  wave  propagating 

through  the  substructure  of  Fig.  1.  In  Fig.  2a  a  rightwardly- 

propagating  stress  wave  of  arbitrary  shape  is  shown  with  its 

wavefront  located  at  cross  section  x,  at  time  t  .  Then,  after  a 

time  duration  of  x  has  elapsed,  the  wave  has  propagated  to  the 

1 2 

right  so  that  its  wavefront  is  located  at  cross  section  x  as 

2 

shown  in  Fig.  2b;  the  shape  of  the  wave  is  unchanged.  After  an 
elapsed  time  of  1  ,  which  is  given  by 


x  -  x 
3  1 


the  wavefront  has  reached  cross  section  x  as  shown  in  Fig.  2c, 

3 

with  the  wave  shape  undistorted  and  the  amplitude  unchanged. 

The  delay  function  related  by  eqn.  (11 )  may  be  given  a  block 
diagram  representation,  as  in  Fig.  3,  where  the  nodes  x  and  x 


reier  to  the  corresponding  substructure  cross  sections,  and 
double-lined  faces  on  the  delay  block  are  used  to  distinguish  it 
from  typical  con-tant-multiplier  blocks.  Block  diagram 
representations  will  be  useful  later  in  visualizing  the  overall 
structure  of  the  control  system  under  study. 

l’ne  operation  represented  by  tne  delay  block  shown  in  Fig.  3 
is  simple;  its  output  signal  or  function  is  merely  the  input 
signal  or  iunction  delayed  by  the  time  constant  with  which  the 
particular  block  is  labelled.  Thus  in  Fig.  3  the  output  stress 

function  c+(x  ,t)  is  merely  the  delay  of  the  input  stress  function 

2 

cr  (x  ,t)  •  Since  the  time  delay  for  wave  propagation  between  x 

i  l 

and  x  is  given  by  eqn .  (12)  as  t  ,  then  eqn.  (11)  describes 

2  12' 

the  operation  of  the  delay  block  in  Fig.  3.  Note  also  that 
eqns .  (12)  and  (13)  may  be  generalized  in  the  following  form: 


(14) 


where  t  is  the  time  delay  for  longitudinal  wave  propagation  in 
eitr.er  direction  between  two  cross  sections  xm  ana  xn . 

Before  addressing  the  control  problem,  two  other  aspects  of 
wave  propagation  in  the  lattice  substructure  will  be  discussed 
and  represented  by  block  diagram  elements,  which  will 
subsequently  be  used  to  Illustrate  the  different  configurations 
and  the  underlying  assumptions  of  a  few  different  controller 
designs . 

First,  waves  may  propagate  through  the  substructure  in 
eitner  the  positive  or  negative  x-directions ,  as  shown  in  eqn. 


11 


I.  — - *T  l  -  T  wr"r  *.  »rv-  «■*■ 


V»  W*  *  .  »  ', 

nv*; 


y.-j* 


(10) . 

Thus,  in  Fig. 4, 

o+(x  ,  c)  is 
i 

shown 

block 

in  the  positive  x 

direction 

( from 

shown 

passing  through 

the  same 

delay 

direction  (from  x  to  x 

). 

The  second  aspect  is  the  principle  of  superposition  in 
linear  systems.  When  two  waves  travelling  in  a  lattice 
substructure  pass  the  same  location,  their  magnitudes  at  that 
cross  section  are  additive.  This  superposition  principle  is 
accountec  for  by  summing  the  values  of  time  functions  for 
leftwardly  and  rightwaraly  propagating  stress  functions  at  a 
cross  section  in  order  to  obtain  the  net  time  function  for  stress 
at  that  cross  section.  Fig.  4  shows  the  standard  block  diagram 
notation  for  such  a  summation  operation.  Just  as  the  net  or 
actual  value  of  stress  at  a  physical  cross  section  is  the  sum  of 
the  instantaneous  values  of  propagating  stress  waves  at  that 
cross  section,  likewise  the  stress  value  at  a  given  cross  section 
in  the  block  diagram  is  the  sum  of  both  the  leftwardly  and 
rightwardly  propagating  stress  functions  for  the  cross  section  at 
that  point  in  time. 


1  2 


X  .  ) 
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1  2 


2  3 


In  Fig.  5>  two  waves  are  shown  approaching  cross  section  x 


a  :a. 
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vV. 
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ve  stress 

■  . 
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for  all 

by  eqn. 

'•*y 

12 


(5).  Fig.  5  shows  both  waves  propagating  at  equal  speed  until 

they  are  fully  superposed  at  time  t  equal  to  t  +  x  .  During 

0  12 

superposition  the  waves  partially  cancel;  their  separate 
magnitudes  are  represented  by  outlines,  and  their  sum  or  net 
magnitude  is  shaded.  After  superposition  the  waves  continue  to 
propagate  in  their  respective  directions  undistorted,  as  Fig.  5 
shows.  Fig.  6  shows  a  block  diagram  model  for  this  length  of  the 
substructure  from  section  x  to  x  . 

l 


a 


III.  FORMULATION  OF  FEEDFOREWARD  CONTROL  MODEL 


A.  Idealized  Controller  Components 

The  principle  of  waveform  superposition  underlies  the  scheme 
to  be  analyzed  xor  vibration  control  in  this  study — namely, 
cancellation  of  propagating  stress  waves  by  generating  and 

superposing  the  negative  of  waveform  disturbances  contained  in 

tne  structure.  The  analysis  will  consider  disturbance  (input) 
waves  which  are  propagating  through  a  lattice  substructure  in  the 
positive  x  direction,  as  well  as  the  operations  of  a  controller 
designed  to  protect  the  region  of  the  substructure  to  the  right 

of  section  x  (refer  to  Fig.  1). 

2 

The  operations  of  the  ideal  feedforeward  controller  to  be 
studied  may  be  described  as  follows.  First,  a  sensor  at  section 
xi  is  used  to  measure  the  amplitude  of  the  stress  wave 
disturbance  as  it  passes  the  section.  The  signal  generated  by 

the  sensor  is  delayed,  and  provides  input  to  the  actuator  at 

cross  section  x  ,  where  the  negative  of  the  incoming  stress  wave 

2 

is  actuated  after  an  appropriate  time  delay.  As  indicated  in 
eqn.  (12),  the  length  of  the  time  delay  will  be  determined  by 

the  distance  i  along  the  substructure  between  the  sensor  and 

1  2 

the  actuator. 

It  is  important  to  clearly  state  the  assumptions  which  are 
involved  in  tne  simplified  modelling  of  the  sensor,  delayer,  and 
actuator  elements  in  this  idealized  controller  design.  They  are 


(l)  Sensor —  The  sensor  is  assumed  to  be  a  transducer  of 
negligible  size  and  mass  which  is  fixed  to  the  substructure 


at  section  x  ,  and  which  generates  a  signal  s(t)  (with  units 
1 

of  volts)  proportional  to  the  longitudinal  stress  at  the 
section  ,  a ( x  , t ) . 

i 

s(t)  *  ~  o(x  ,t)  (15) 

Gi  1 

here  (1/G  )  is  the  constant  of  proportionality  (or  "gain") 
i 

associated  with  tne  sensor,  and  has  units  of  (volt/psi). 
The  actual  operating  principle  (piezoelectric, 

electromechanical,  strain- resistive ,  etc.)  is  not  considered 
in  the  present  analysis. 

(2)  Delayer —  The  signal  s(t)  is  fed  to  a  signal  delay  element, 
whose  output  s'(t)  is  given  by 


s'(t)  =  s(t-r'  )  ;  for  c  >  r'  (16) 

12  —  1 2 


where 

T  '  is 

12 

the  time  delay 

introduced  by 

the 

delayer. 

In 

the 

case  of 

ideal  operation 

,  x  '  will  equal 

1  2 

x  ,  which 
1  2 

is 

the 

actual 

transportation 

lag  of  stress 

wave 

s  propagat 

ir  - 

■*.  o 

frcr;. 

section 

x  to  section  x 

• 
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(3)  actuator—  The  actuator  is  located  at  a  particular  section 

in  the  rod,  x  .  It  is  approximated  by  a  disk  whose  thickness 

2 

is  negligible  and  whose  impedance  is  matched  to  that  of  the 
rod,  so  that  stress  waves  propagate  past  cross  section  x 

2 

unaffected  by  the  presence  of  the  actuator. 

It  is  assumed  that  the  controller  can  apply  a  time  varying 
axial  force  F(t)  to  the  disk,  and  that  the  value  of  F(t)  is 


proportional  to  the  signal  s'(t)  received  from  the  delayer. 

F(t)  =  G  •  s '  ( t )  (17) 

2 

Here,  G2  is  the  actuator  gain,  and  has  units  of  (lb/volt). 

Cnee  again,  the  present  analysis  does  not  consider 
whether  the  force  is  actually  generated  piezoelectrically , 
e lec tromagne tica 1 1 y ,  etc.  What  is  important  is  that  the 
controller  is  assumed  to  generate  time-varying  uniform  stress 
distributions  across  the  faces  of  the  substructure. 

b.  Description  of  Actuator  Operation 

Fig.  7a  shows  the  idealized  actuator-substructure  assembly 
intact,  where  the  rod  has  cross  sectional  area  A.  In  Fig.  7b 
segments  of  the  lattice  substructure  have  been  isolated  as  free 
boay  diagrams,  to  show  the  applied  forces  and  resultant  stresses. 
In  Figs.  7a  and  7b,  the  positive  x  direction  is  indicated  and 
tensile  stresses  are  defined  positively  as  shown. 

It  is  necessary  to  derive  expressions  for  the  applied  stress 
functions  0(t)  and  ^(t)  in  the  substructural  element.  Examining 
Fig.  7b,  a  balance  of  forces  can  be  written  for  the  massless 
actuator  component. 

ZF  =  0  =  F(t)  +  A-o  (t)  -  A*a  (t)  (18) 

X  2  1 

Therefore,  F(t)  =  A[ai  (t)  -  Oz  (t)  ]  (19) 

Geometric  compatibility  requires  that  displacements  and 
velocities  at  the  actuator  interface  be  equal.  Thus, 


u  (t)  *  u  (t) 
1  2 


3u  (t)  3u  (c) 
1  _  2 _ 


Iiext,  the  cross  sectional  velocities  can  be  expressed  in  terms  of 
the  stresses  on  the  corresponding  faces  [2]. 

3u  (t)  c 

— -  =  —  o  (t)  (22) 

3t  E  1 

3u  (t)  c 

-  =  -Z-O(t)  (23) 

3t  E  2 

Combining  eqns.  (20)  through  (23)  yields  the  requirement  on 
the  stresses  as 


—  a  (t) 
E  i 


—  a  (t) 
v  ? 


Therefore, 


a  (t)  =  -a  (t) 
1  2 
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This  result,  together  with  the  force  balance  relation,  (eqn. 
v19))»  gives  the  expressions  for  the  resultant  stress  functions 
at  me  actuator  in  terms  of  the  forcing  function  as 


a  (t)  = 
x 


o  (t) 
2 


t  •  -  *  W  " 


Finally,  it  is  necessary  to  write  expressions  for  the 
propagation  characteristics  of  the  stress  waves  which  the 
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From  the  end-actuation  analogy,  the  stress  waves  generated 
on  the  substructure  face  to  the  right  cf  the  actuator,  o (t), 
must  propagate  into  this  "semi-infinite  half"  of  the 
substructure;  that  is,  they  must  propagate  in  the  positive  x 
direction.  Thus,  the  equation  for  the  propagation  of  controller- 
induced  stress  waves  in  the  portion  of  the  rod  to  the  right  of 
the  actuator  can  be  written  by  selecting  the  appropriate  term 
from  tne  D'Alembert  solution  to  the  wave  equation. 


a+(x,t)  = 


for  x  >  x  ,  t  >  0  (28) 

2  ~ 


The  superscript  "  +  " 
waves  in  the  positive  x 
function  F  in  eqn.  (28) 
function  continues  to 
direction  at  a  speed  c 


denotes  the  propagation  of  these  stress 
direction.  Note  that  the  argument  of  the 
is  in  units  of  time,  not  position;  the 
represent  a  wave  travelling  in  the  x 


o 

Likewise,  actuated  stress  waves  which  are  generated  on  the 
substructure  face  to  the  left  of  the  actuator  must  travel  to  the 
left,  again  based  on  the  division  of  the  problem  into  two  cases 
of  end-actuation  of  semi-infinite  substructures.  Thus,  the 
equation  for  these  lef twaroly-propagating  stress  waves  also  takes 
the  form  of  the  appropriate  term  from  the  D'Alembert  solution  as 
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cT(x,t)  =  —  •  F  (  - +  t  ]  ;  for  X  <  X  ,  t  >  0  (29) 

2A  \  Co  /  2  ” 


Together,  eqns.  (28)  and  (29)  summarize  the  results  for 
propagating  stress  waves  which  are  generated  by  the  actuator,  and 
they  express  the  waves  in  terms  of  the  value  of  the  actuator 
forcing  function,  F(t).  Using  eqns.  (15)  to  (17),  eqns.  (28)  and 
(2S)  for  the  controller-inducea  stress  waves  may  be  expressed  as 
functions  of  the  stress  detected  by  the  sensor  at  section  x 


a+(x,t) 


for  x  >  x  ,  c  >  t' 

2  —  12 


(30) 
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for  x  <  x  ,  t  >  x' 

2  —  1 2 

tig.  8  presents  a  schematic  of  the  actuator  output  response  to 
the  arbitrary  actuator  input  signal  s'(t)  snown.  Fig.  9  shows  a 
modified  blocx  diagram,  composed  of  the  elements  discussed  in 
Section  II,  which  models  the  general  response  of  the  controlled 
substructure  to  a  stress  wave  disturbance  originating  at  xq ,  with 
the  controller  formulated  as  above.  In  Fig.  9*  solid  arrows 

continue  to  represent  the  "flow"  or  travel  of  stress  wave 
functions,  while  dashed  lines  have  been  used  to  represent 


(electrical)  signal  flow  in  the  controller 


IV.  RESPONSE  OF  THE  CONTROLLED  SUBSTRUCTURE 


It  is  now  possible  to  analyze  the  response  of  the  controlled 
lattice  substructure  to  an  arbitrary  input  disturbance  at  x 

o 

Two  important  features  characterize  the  response  of  the  system 
as  it  is  presently  formulated. 

The  first  nay  oe  likened  to  the  "stress  intensification" 
pnenomenon  which  is  characteristic  of  stress  wave  reflection  frc 
s  fixed  boundary  (such  as  a  rigid  wall).  As  part  of  a  wave 
reflects  back  from  a  fixed  boundary,  it  may  be  superimposed  with 
the  still-incoming  portion  of  the  wave,  so  that  higher  amplitude 
transient  waveforms  result.  The  same  effect  iB  the  result  of 
the  superposition  of  incoming  waves  approaching  the  actuator  at 

x  from  x  and  the  already-actuated,  lef twardly-propagating 

10 

stress  waves,  which  are  duplicates  of  previously  incoming  waves. 

The  series  of  pulse  diagrams  for  a(x,  t)  in  Fig.  10 
illustrate  this  phenomenon.  In  Fig  10,  stress  pulses  with  dashed 
outlines  represent  stress  waves  generated  by  the  controller 
actuator.  Pulses  with  solid  outlines  are  those  of  incoming 
stress  wave  disturbances.  Consistent  with  all  stress  pulses 
presented  in  this  study,  the  shaded  portion  of  a  pulse  represents 
the  net,  or  physically  realized  magnitude  of  the  stress  waves  in 
the  substructure. 

Thus,  in  Figs.  10c  through  lOf,  the  incoming  stress  wave 
(solid  outline)  is  completely  cancelled  by  the  rightwardly 
propagating  controller-generated  stress  wave  (dashed  outline,  for 
x  >  x2  ),  so  that  both  are  left  unshaded,  indicating  that  their 
net  magnitude  is  zero.  They  superpose  and  cancel,  so  that  the 


portion  01  the  lattice  substructure  to  the  right  of  the 

controller  is  isolated  from  any  stress  disturbance  fcy  the 

incoming  wave.  Note  that  in  Fig.  10  T  J  is  assumed  to  be  equal 

12 

to  T 

1  2 

In  Fig.  10c  the  phenomenon  of  stress  intensification  is 
shown.  The  lef twardly-propagating  stress  wave  generated  by  the 
controller  is  momentarily  superposed  on  a  portion  of  the 

incoming  wave  disturbance  which  is  travelling  to  the  right  past 
section  ;  the  result  of  this  superposition  is  transient  stress 
values  which  are  of  higher  magnitude  than  either  cf  the  two 

waves’  peak  values  alone. 

There  is  a  second  problematic  feature  of  the  controller. 
Note  that  in  Fig.  10ef  the  lef twardly-propagating  stress  wave 
(which  is  generated  by  the  actuator  as  it  also  produces  the 

negative  of  the  incoming  wave),  is  passing  the  sensor  located  at 
cross  section  X1  •  As  the  controller  is  presently  formulated, 
both  incoming  ( rightwardly-propagating )  stress  waves  and 
outgoing  ( lef twardly-propagating )  stress  waves  are  sensed  by  the 
controller  as  "substructure  disturbances".  That  is,  the  sensor 
currently  generates  a  signal  (delayed  for  input  to  the  actuator) 
based  on  tne  instantaneous  magnitude  of  any  stress  wave  passing 
tne  cross  section  x^  ,  regardless  of  the  wave's  source. 

For  example,  the  stress  waves  shown  in  Fig.  lOg  are  being 
generated  by  the  actuator  as  a  result  of  the  outgoing  stress  wave 
which  is  shown  passing  the  sensor  earlier  in  Fig.  lOe.  Note  that 
the  leftward-propagating  wave  generated  in  Fig.  lOg  will  in  turn 
pass  the  sensor.  Thus,  a  single,  non-repeated  incoming  stress 


pulse,  though  cancelled  perfectly  by  the  controller,  will 
nonetheless  initiate  a  process  of  endless  controller  "echoing", 
as  lei twaraly-propagating  controller-generated  stress  waves  are 
repeatedly  sensed,  producing  new  actuations. 

These  two  problematic  features  of  the  controller,  namely 
stress  intensification  (Fig.  10c)  and  "echoing",  (Figs.  9g  and 
9n)  can  be  compounded  to  the  degree  that  the  system  becomes 
unstable  in  response  to  disturbances. 

Fig.  11  shows  how  such  instability  can  arise  in  response  to 
a  simple  step  stress  disturbance.  Cashed  lines  again  represent 
the  outlines  of  controller-generated  stress  waves.  The  dashed 
outline  cf  these  actuated  waves  is  shown  slightly  below  the  solid 
outline  of  the  incoming  stress  wave  for  clarity  only;  such  a 


representation  is 

intended 

to 

indicate 

equal  magnitude 

nonetneless.  As  in 

Fig.  10, 

Fig. 

11  assumes 

that  t  '  =  t 

12  12 


Fig  11b  shows  that  'stress  intensification",  or  additive 
superposition,  has  begun  to  occur  by  t  equal  to  t  ♦  Ti z/2,  so  that 
a  wave  of  twice  the  amplitude  of  the  incoming  wave  results.  As 
tne  outgoing  actuated  wave  propagates  further  to  the  left,  so 
does  the  additive  effect  of  its  superposition  with  the  incoming 
wave.  By  Fig.  lid,  when  t  equals  tfl  +  (3/2)f.i2  ,  the  higher 
amplitude  wave  has  passed  the  sensor  at  cross  section  . 
Accordingly,  as  snown  in  Fig.  Ilf,  the  actuator  begins  to 
generate  a  stress  wave  of  twice  the  incoming  wave's  amplitude  at 
a  time  which  is  3TJ2after  the  incoming  wavefront  initially  passed 
the  sensor. 

Figs.  Ilf  and  11g  show  that  as  the  generated  waves  continue 
to  superpose  with  the  incoming  wave,  then  cancellation  is  no 


longer  achieved  to  the  right  of  x  ,  because  the  controller 

2 

generates  waves  which  are  twice  the  necessary  amplitude. 
^  Finally,  Fig.  11g  further  shows  that  higher  amplitude  waves  will 

continue  to  pass  the  sensor,  leading  to  generated  waves  of  still 
greater  amplitude.  The  system  is  unstable  in  response  to  a  step 
>  input. 

Fig.  12  illustrates  the  case  for  a  step  sinusoidal  input. 
Again,  the  lef twaraly-propagating  waves  which  are  generated  by 
|  the  controller  las  it  operates  to  cancel  the  incoming  waves) 

superpose  with  the  incoming  stress  waves.  Stresses  of  double  the 
peak  amplitude  of  the  incoming  wave  are  subsequently  measured  by 

'  the  sensor  at  x  (see  Fig.  12e),  and  the  controller  is  required 

1 

to  generate  waves  of  increasingly  higher  amplitude  (dashed 
outlines  in  Figs.  12g  and  12h) .  Thus  the  system  is  potentially 
&  unstable  in  response  to  a  step  sinusoidal  input. 


vs 


V.  COiiTkOLLEK  STABILIZATION 


A  preliminary  approach  lor  actively  controlling  the 
propagation  of  longitudinal  vibrations  in  an  idealized  lattice 
substructure  has  been  discussed  in  Sections  II  through  IV.  As 
developed  in  Section  IV,  the  controlled  substructure  system  is 
unstable  in  response  to  input  disturbances.  The  instability 
arises  because  the  sensor  responds  to  both  incoming  and 
controller-generated  lef twardly-propagating  stress  waves 

indiscriminately.  It  is  desirable  to  modify  the  controller  to 
prevent  tnis  instability. 

One  solution  is  to  redefine  the  operational  characteristics 
of  the  sensor.  It  is  reasonable  to  assume  a  sensor  which  is 
sensitive  only  to  waves  propagating  in  the  desired  direction — in 
this  case,  the  positive  x  direction.  Employing  such  a  sensor, 
the  controller  would  not  be  reactivated  by  the  leftwardly- 
propagating  stress  waves  which  the  actuator  generates,  and  thus 
the  system  would  not  exnibit  the  instability  characteristics 
discussed  in  Section  IV. 

A  second,  more  general  solution  to  the  problem  of 
instability  would  be  to  add  a  loop  to  the  controller  which 
effectively  deducts  the  delay  of  the  actuator  output  from  the 
sensor's  output.  Implementation  of  this  solution  is  illustrated 
by  a  block  diagram  representation  in  Fig.  13.  As  shewn  in  Fig. 
13,  the  signal  which  provides  input  to  the  actuator  is  also  fed 
back  to  be  subtracted  from  the  sensor  output,  through  a  delayer 

labeled  x  ' ' .  Since  the  objective  of  this  additional  loop  is 

1 2 

to  deduct  from  s(t)  that  portion  of  the  signal  which  corresponds 

24 


to  the  leftwarcly-propagating  (actuator-generated)  portion  of 

0 ( x  »t),  then  it  is  necessary  that  t  ’  '  be  exactly  equal  to  x 

1  12  12 

in  order  to  totally  cancel  the  actuator-generated  portion  of 

a ( x  ,t).  A  nonzero  difference  between  x  "  and  x  is  thus  an 

1  12  12 

additional  possible  source  of  controller  error;  the  effects  of 
such  a  nonzero  difference  will  not  be  analyzed  here,  but  the 
study  of  such  effects  constitutes  an  appropriate  topic  fcr 
further  research. 

For  the  remaining  analysis,  the  successful  application  of 
this  latter  more  general  stabilized  method  to  the  controller  will 
be  assumed.  Thus,  in  the  remaining  analysis,  it  is  assumed  that 

x  ' '  equals  x  ,  so  that  the  deduction  of  the  delayed  value  of 

12  12 

s'(t)  irom  the  sensor  output  occurs  at  precisely  the  same  instant 
that  the  sensor  is  measuring  the  lef  twardly-propagatir.g 
disturbance  which  was  generated  in  response  to  s'(t)  and  which 

nas  travelled  back  tnrough  the  substructure  from  x  to  x  .  The 

2  1 

response  of  the  controlled  lattice  substructure  to  a  step  input 
disturbance,  following  this  controller  modification,  is  shewn  in 
Fig.  14.  In  Fig.  14,  the  effect  of  stress  intensification  has 

propagated  back  past  the  sensor  at  section  x  by  time  t  equal  to 

1 

t  +  2x  .  However,  as  Fig.  14  shows,  for  time  t  equals  3  X'  and 
beyond,  the  actuator  never  doubles  the  value  of  the  controller- 
generated  wave,  because  a  value  associated  with  the  leftwardly- 
propagating  disturbance  is  effectively  subtracted  from  the  sensor 
output,  as  schematized  in  Fig.  13.  Fig.  14  may  be  compared  with 
Fig.  11,  discussed  in  Section  IV,  which  illustrates  the  response 
of  tr.e  unmodified  controller-substructure  assembly  to  an 
identical  step  input  cisturoance. 


VI 


INPUT/OUTPUT  CHARACTERISTICS  OF  STAELIZED  FEEEFOREWARD- 
CONTRCLLED  SUBSTRUCTURE 


A.  Transfer  Function  Magnitude 


It  is  useful  to  characterize  a  linear  system  by  its  transfer 
function,  or  (steady  state)  frequency  response,  H(w),  which  is 
given  by 


H(w) 


Y(u) 

X(u) 


(32) 


wnere  *>  denotes  raaian  frequency,  X(ui)  is  a  steady  harmonic 
input  to  tne  system  and  Y(u)  is  the  system  output.  The  transfer 
function  for  the  feecf oreward-ccntrolled  substructure  which  has 
been  treated  in  the  previous  sections  is  derived  in  Appendix  E, 
and  found  to  be  given  by 


Hf  (w) 


-iOJT 

e 


1  2 


Ke-M>’ 


(33) 


where  x  and  T  ’  are  the  time  delays  for  wave  propagation 
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through  the  substructure  and  the  controller,  respectively;  the 
controller  net  feedforeward  amplification,  given  by  Gz/2AG1  in 
Section  III,  has  been  replaced  by  K  for  simplicity;  and  i  =  * -1  . 

It  can  be  seen  from  eqn.  (33)  that  the  transfer  function  is 
a  complex  function.  Defining 


£ 


-  T 

1  2 


(34) 


where  e  denotes  the  difference  between  the  time  delay  of  waves 

propagating  in  the  substructure  x  and  the  delay  of  the 

1  2 

controller  t  '  ,  then  Appendix  B  shows  that  the  square  of  the 

1  2 

transfer  function  magnitude  is  given  by 


(35) 


2 

( Hf  (ui)  |  =  1  -  2Kcos (cue)  +  K2 

In  terms  of  the  controlled  substructure  system  under  study, 
the  magnitude  of  the  transfer  function  given  in  eqn.  (35) 
represents  the  complex  magnitude  of  the  ratio  of  the  system 
output  to  the  system  input  where  both  are  expressed  as  complex 
functions  in  the  frequency  domain.  In  the  time  domain,  the 
system  input  and  output  are  assumed  to  be  stress  functions 
o+(xi  ,t)  and  °(x3  »t) ,  respectively,  both  are  representable  by 
harmonic  components,  x  and  x  denote  the  locations  of  cross 

1  3 

sections  on  the  substructure  as  shown  in  Fig.  9,  and  the 
superscript  "+"  indicates  that  only  stress  waves  propagating  in 
the  positive  x  direction  past  section  x  i  are  actually  input  into 
the  system. 


As  long  as  the 

time 

domain  input  and 

output  of 

the 

controlled 

substructure 

are 

sinusoidal  stress 

functions 

with 

frequency 

,  a  transfer 

function  magnitude  i 

Kf(<A)l  of 

less 

than  unity 

represents 

system  attenuation  of  an 

input  stress 

function.  A  transfer  function  magnitude  equal  to  zero  represents 
complete  attenuation  or  cancellation  of  the  input  stress  function. 
A  transfer  function  magnitude  |Hf(w)|  greater  than  unity 
indicates  controller  operation  which  effectively  amplifies 
incoming  sinusoidal  stress  functions  of  frequency  w  ,  so  that  the 
amplitude  of  the  stress  wave  output  to  cross  section  x^  in 
Fig.  9  is  greater  than  that  of  the  input  stress  wave. 


b.  Operational  Ranges  of  Controller  Error 

Two  types  of  controller  error  are  considered  in  this  study: 

(1)  nonzero  values  for  the  difference  between  controller  and 
substructure  delays,  e  ,  which  is  given  by  eqn.  (34);  and 

(2)  values  of  the  net  feedforeward  controller  gain  K  not  equal  to 
the  substructure  gain  of  unity. 

There  are  also  two  possible  operational  goals  for  the 
feedforeward  controller,  within  the  broad  domain  of  "active 
vibration  control".  The  first  possible  goal  is  "ideal” 
controxler  performance— achieving  complete  superpositions  1 
cancellation  of  the  incoming  stress  wave  disturbance ,  thus 
completely  isolating  cross  section  from  any  ( rightwardly- 

propagating)  vibrational  disturbances  (see  Fig.  9).  The  second, 
more  modest  goal  of  controller  implementation  would  be 
achievement  of  disturbance  attenuation,  if  not  full  cancellation. 
In  this  case  the  objective  might  be  an  output  stress  wave  whose 
amplitude  is  some  fraction  of  the  input  disturbance  amplitude, 
thus  requiring  a  transfer  function  magnitude  less  than  unity. 
The  present  subsection  identifies  the  allowable  ranges  of 
controller  error  for  achievement  of  each  of  these  two  controller 
objectives . 

(1)  Achieving  complete  disturbance  cancellation 

"Ideal"  operation  of  the  controller  calls  for  complete 
cancellation  of  incoming  stress  wave  disturbances  which  are  input 
to  the  controlled  substructure  system.  It  is  possible  to  derive 


the  requirements  on  the  controller  delay  error,  €  ,  given  by  eqn. 
(34),  and  the  requirements  on  the  net  feedforward  controller 
gain  K. 

Eqn.  (35)  gives  the  expression  for  the  square  of  the 
system  transfer  function  magnitude.  Achievement  of  full  input 
disturbance  cancellation  requires  a  system  output  of  zero,  and 
thus  a  transfer  function  magnitude  of  zero.  Kith  this 
requirement  on  |  H  ( to)  |  ,  eqn.  (35)  becomes 

K2  -  2Kcos(we)  +1  =  0  (36) 


Eqn . 
by 


(36)  is  a  quadratic  equation,  so  that  its  solution  is  given 


2cos(toe)  ± 


d 


■  2cos  (<oe) 


-  4 


(37) 


which  may  be  simplified  and  rewritten  as 


K 


cos  (toe) 


cos2 (we)  -  1 


(38) 


Eqn , 

(38) 

is 

the 

requirement 

for 

full 

controller 

cancellation 

of 

incoming 

disturbances  . 

If 

K ,  the 

controller 

feedforward 

gain 

,  is 

to 

take  on  only 

real 

values , 

then  the 

difference  under  the  radical  in  eqn.  (38)  must  be  nonnegative. 
This  requirement  states  that 


cos2 (we)  >  1 


(39) 


ana  tnus  ideal  controller  operation  requires  that 


toe  *  irn  (40) 

where  n  is  any  integer,  n  =  ...-2, -1,0,1 ,2.... 

Using  the  relation 

2tt 

T  =  —  (41) 

to 

where  1  is  the  period  (in  seconds)  of  a  sinusoidal  waveform  with 
radian  frequency  to  ,  the  requirement  on  the  controller  delay 
error  (eqn.  (40))  can  be  rewritten 


e  =  j  (nT) 


Eqn.  (42)  indicates  that  for  complete  cancellation  of  a 
sinusoidal  input  disturbance,  the  duration  of  the  controller 
delay  error  e  must  be  equal  to  some  integer  multiple  of  half  the 
input  waveform  period  T. 

Next,  combining  eqns.  (40)  and  (38)  yields 

!1  for  n  even 

(43) 

-1  for  n  odd 

which  is  the  second  requirement  on  the  controller  parameters 
which  must  be  fulfilled  if  the  controller  is  to  operate  ideally. 

Thus,  in  summary,  for  the  controller  to  operate  such  that 
complete  cancellation  of  input  disturbances  is  achieved,  both 
eqns.  (42)  and  (43)  must  be  satisfied  simultaneously  by  the 
controller  parameters. 
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(2)  Achieving  disturbance  attenuation 


In  the  case  where  the  goal  of  controller  implementation  is 
achieving  attenuation  of  stress  wave  disturbances  which  are  input 
to  the  system,  it  is  necessary  that  for  an  input  sinusoidal 
stress  function  of  radian  frequency  oj  ,  the  magnitude  j  H f  (uo )  |  of 
the  system  transfer  function  must  be  less  than  unity.  Solving 
eqn.  (35)  for  this  requirement  yields 

1  >  1  -  2Kcos (we)  +  K2  (44) 


Eqn.  (44)  may  be  simplified  to  obtain 


K  <  2cos(oje)  ;  for  K>0 

K  >  2cos(o)e)  ;  for  K<0 


(45) 


Eqn.  (45)  is  plotted  in  Fig.  15»  where  all  value  pairs  for  K 
and  cue  which  satisfy  the  inequality  lie  in  the  shaded  portion  of 
tne  plane,  but  do  not  include  the  points  where  the  line  K=2cos(coe) 
crosses  the  (uejaxis. 

Examination  of  Fig.  15  and  eqn.  (45)  indicates  that  while 
there  are  no  upper  or  lower  bounds  for  acceptable  values  of  the 
parameter  (uje),  the  magnitude  of  the  controller  feedforward  gain 
K  must  never  reach  or  exceed  a  value  of  2  if  attenuation  is  to 
be  achieved.  Also,  the  parameter  (ue)  must  not  equal  any 
odd  integer  multiple  of  it/ 2,  and  the  parameters  K  and  (toe)  cannot 
be  set  independently,  if  attenuation  is  to  be  achieved. 
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C.  tnysical  Basis  for  Frequency  Dependence  of  Transfer  Function 
Magnitude 

£qn.  (35) »  the  expression  for  the  magnitude  of  the 

controlled  substucture  system  ;  ansfer  function,  is  plotted  in 

2 

Fig.  E-3  as  the  transfer  function  magnitude  squared,  |Hf(w)|  , 

versus  the  parameter  we  .  Eqn.  (35)  and  Fig.  B-3  show  that  the 
magnitude  of  the  transfer  function  depends  on  the  frequency  w  of 
the  incoming  disturbance.  That  is,  for  a  constant  feedforward 
gain  K  and  a  constant  difference  e  between  substructure  and 
controller  delays,  the  magnitude  of  the  system  transfer  function 
will  be  different  for  different  input  signal  frequencies. 
Therefore,  since  actual  inputs  to  the  controlled-substructure 
system  may  be  waveforms  composed  of  several  harmonics  at 
different  frequencies,  the  magnitude  of  certain  such  harmonics 
will  be  enhanced  in  the  system  response  while  that  of  harmonics 
at  other  frequencies  will  be  diminished.  In  fact,  as  Fig.  B-3 
illustrates,  the  magnitude  of  the  transfer  function  is  periodic 
in  we  .  That  is,  as  we  is  increased  (or  decreased)  by  2u  , 
the  magnitude  of  the  transfer  function  repeats  itself. 

To  understand  physically  the  dependence  of  the  system 
transfer  function  magnitude  on  the  input  waveform's  harmonic 
frequencies,  it  is  useful  to  introduce  the  concept  of  a  waveform 
period.  Making  use  of  eqn.  (41 )»  where  T  is  the  period  (in 
seconds)  of  a  waveform  whose  radian  frequency  is  w  ,  note  that 

we  =  (46) 
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Substituting  eqn.  (46)  into  eqn.  (35)  then  gives 


[Hf(2ir/T)j  =  1  -  2Kcos 


[»«)] 


Eqn.  (47)  is  plotted  in  Fig.  16  to  illustrate  that  the  magnitude 
of  the  transfer  function  is  periodic  in  (e/T).  It  is  instructive 
to  consider  That  the  ratio  (e/T)  is  a  measure  of  e  as  a 
percentage  of  the  waveform  period.  For  a  waveform  with  constant 
period  T,  Fig.  16  shows  that  as  the  difference  in  the  controller 
and  substructure  time  delays  e  is  increased  from  0  to  50  percent 
of  T,  the  magnitude  of  the  transfer  function  is  increased  from 
its  minimum  to  its  maximum  value. 

As  developed  in  Section  VI-B,  only  when  the  magnitude  of  the 
net  controller  gain  equals  unity  is  it  possible  for  the 
controlled  system  to  completely  attenuate  incoming  signals. 
Figs.  16  and  E-3  present  the  frequency  dependency  of  the  transfer 
function  magnitude  for  the  general  case  where  K  is  unspecified. 
The  remaining  analysis  will  assume  K  equal  to  unity,  (zero  gain 
error) ,  and  will  quantify  the  effects  of  controller  delay  error 
given  this  assumption. 

Fig.  17  is  a  plot  of  the  square  of  the  system  transfer 
function  magnitude  versus  the  parameter  (e/T),  with  K  set  equal  to 
unity.  Of  particular  interest  is  the  dashed  horizontal  line 
in  Fig.  17  representing  a  transfer  function  magnitude  of  unity. 
Above  tnis  line  the  controlled  system  amplifies  an  input 

disturbance;  below  it  the  input  is  attenuated. 

Eqn.  (45)  may  be  solved  for  K  equal  to  unity  to  determine 
the  ranges  of  the  parameter  (  e  / T )  over  which  the  system 
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attenuates  an  input.  Eqn.  (45)  becomes 

cos(we)  >  1/2  (48) 

which  can  be  rewritten,  using  eqn. (41),  as 

cos  >  1/2  (49) 

Eqn.  (4S)  is  solved  for  (e/t)  in  the  following  ranges: 

(n  -  1/6)  <  -  <  (n  +  1/6)  (50) 

T 

wnere  n  is  any  integer. 

Eqn.  (50)  thus  specifies  tl  i  ranges  of  the  parameter  (e/T) 
for  which  the  feedforward  controlled  substructure  system 
successfully  attenuates  an  incoming  stress  wave  disturbance  for  K 
equal  to  unity. 

When  the  feedforward  gain  is  held  constant  at  unity,  then  for 
(e/T)  in  the  ranges  (n-1/6  to  n+1/6)  where  n  is  any  integer,  the 
system  attenuates  a  steady  sinusoidal  input  signal.  Likewise, 
for  (e/T)  in  the  ranges  (n+1/6  to  n+5/6)  the  system  amplifies  a 
steady  sinusoidal  input  signal. 

tiathmatically ,  this  result  is  found  by  solving  eqn.  (46) 
when  K  and  I H ^ ( 2  VT)|  are  both  equal  to  unity.  Physically,  it 
arises  because  for  ( e/t)  in  the  ranges  (n-1/6  to  n+1/6)  the 
sinusoidal  output  of  the  controller  and  substructure  tend  to 
subtract  in  phase,  while  for  (e/T)  in  the  ranges  (n+l/6  to  n+5/6) 
the  controller  and  substructure  output  tend  to  add  in  phase. 

Fig.  18  presents  graphically  the  relationship  discussed 
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above  for  £  as  a  varying  percentage  of  a  constant  period  T,  and 
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K  equal  to  unity.  Fig.  18  shows  the  stress  at  location  x  as  a 
function  of  time  for  a  sinusoidal  input  of  unit  amplitude.  In 
Fig.  18a  e  =  0  ,  so  that  the  substructure  output  stress 

function  at  cross  section  x  (solid  line)  and  the  controller 

2 

output  stress  function  at  cross  section  x  (dashed  line)  subtract 

2 

completely  (that  is,  they  are  exactly  160  degrees  out  of  phase); 
ana  the  magnitude  of  their  sum  is  equal  to  zero  everywhere.  This 
is  the  case  of  ideal  controller  operation  (feedforward  gain  K 
equal  to  unity  and  the  controller  delay  exactly  equal  to  the 
substructure  delay),  in  which  complete  superpositions  1 
cancellation  is  achieved. 

In  Fig.  18b  (e/T)  equals  1/6,  and  the  system  output  waveform 
(the  sum  of  the  controller  and  substructure  outputs)  is  indicated 
by  the  shading.  Fig.  18b  shows  that  the  peak  amplitude  of  the 
system  output  is  equal  to  the  peak  amplitude  of  the  substructure 
output  (and  thus  that  of  the  input  signal,  since  the  substructure 
gain  is  equal  to  unity),  wnich  indicates  that  the  magnitude  cf 
the  transfer  function  is  unity  for  K  equal  to  unity  and  for  (e/T) 
equal  to  1/6. 

In  Fig.  18c,  (e/T)  equals  1/2.  Substructure  and  controller 
output  add  in  phase  in  Fig.  16c,  yielding  an  output  waveform 
whose  peak  amplitude  is  twice  that  of  the  substructure  output, 
and  thus  is  also  twice  that  of  the  input  waveform.  In  this 


case  then,  the  magnitude  of  the  transfer  function  is  equal  to  2, 
which  is  identical  to  the  result  shown  for  k'  equal  to  unity  and 
(e/T)  equal  to  1/2  in  Fig.  17  and  eqn.  (47). 


Finally,  Fig.  18d  shows  the  system  response  when  (  e/T)  equals 
unity.  The  results  are  indistinguishable  from  those  in  Fig.  18a, 
illustrating  the  periodicity  of  |H  (2ir/T)|  with  respect  to  (e/T). 
lihen  the  controller  output  waveform  leads  or  lags  the 
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substructure  output  waveform  by  e  equal  to  any  integer  multiple 
of  the  period  (in  this  case  the  integer  is  unity)  then  the  two 
output  waveforms  are  superposed  exactly  18C  degrees  out  of  phase, 
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and  thus  they  completely  cancel. 

In  actual  operation  of  the  controlled  substructure  system, 
it  is  more  likely  that  the  delay  error  e,  given  by  eqn.  (34), 
will  be  fixed  and  the  system  will  be  subject  to  an  input  composed 
of  numerous  narmonics  which  are  distributed  over  a  range  of 
frequencies.  The  magnitude  of  the  transfer  function  for  each 
harmonic  will  vary  according  to  the  particular  harmonic's 
frequency  in  a  manner  dictated  by  eqn.  (47)  and  identical  to  that 
discussed  above  ana  presented  in  Fig.  18.  The  parameter  l/T 
varies  linearly  with  the  radial  frequency  u>  ,  as  eqn.  (41) 


l  3 


indicates,  so  that  for  higher  frequencies  the  period  T  is 
shortenea  ana  a  constant  e  becomes  a  larger  percentage  of  T, 


increasing  the  ratio  (e/T). 


The  consequent  frequency 


aepender.ee  of  the  transfer  function  magnitude  is  presented 
grapnically  in  Fig.  19,  for  which  e  is  constant,  K  is  held  equal 
to  unity,  ana  the  input  waveform's  frequency  (and  thus  its 
period)  are  varied. 

When  a  waveform's  period  is  equal  to  6 e  ,  then  (e/T)  equals 
1/6,  and  |  H  (2  tt/T)|  is  equal  to  unity  (for  K  equal  to  unity). 
This  case  is  shown  in  Fig.  19a,  where  the  substructure  and 
controller  outputs  are  represented  by  solid  and  dashed  lines, 
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respectively,  and  their  sum  (the  system  output)  is  shaded 


In  Fig.  19b*  the  frequency  has  increased  to  3  times  its 
value  in  Fig.  19a,  reducing  the  period  to  twice  the  length  of  e, 
so  that  the  controller  and  substructure  outputs  are  exactly  in 
phase.  In  this  case  the  controller  and  substructure  outputs 
additively  superpose  to  produce  a  system  output  whose  magnitude 
is  twice  the  magnitude  of  the  input  (that  is,  twice  the 
magnitude  cf  the  unity-gain  substructure's  output).  This 
transfer  function  magnitude  of  i  H f  ( 2  tt  / T )|  equal  to  2  is  exactly 
the  magnitude  specified  in  Fig.  17  for  (e/T)  equal  to  1/2  and  K 
equal  to  unity. 

Finally,  in  Fig.  19c,  the  frequency  is  6  times  its  value  in 
Fig.  19a,  so  that  T  now  equals  £  exactly,  and  the  controller  and 
substructure  outputs  are  exactly  half  a  period  out  of  phase  so 
that  they  cancel  entirely.  This  case  illustrates  complete  input 
attenuation,  (|Hf(27T/T)|  equals  0)  ,  which  is  what  eqn.  (47)  and 
Fig.  17  give  when  K  equals  unity  and  ( e/ T )  takes  on  an  integer 
value . 

Together,  Figs.  18  and  19  provide  several  system  output 
examples  which  illustrate  the  physical  basis  for  the  frequency 
dependence  of  the  controlled  substructure  transfer  function 
magnitude.  This  dependency  is  derived  in  Appendix  B,  and  is 


presented  in  eqns.  (35)  and  (47),  and  in  Figs.  E-3  and  16 


VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


To  avoid  instability,  a  feedforward  controller  used  to 
actively  control  wave  propagation  in  lattice  substructures  should 
be  designed  to  effectively  sense  and  cancel  only  incoming 
disturbances.  This  may  be  achieved  by  utilizing  a 
unidirectional  sensor;  a  more  generally  applicable  solution  is  to 
ada  a  correction  feedback  loop  in  the  controller. 

Two  possibilities  for  controller  error  include:  (l)  a  net 
feedforward  controller  gain  that  does  not  match  that  of  the 
substructure;  and  (2)  a  controller  delay  time  which  differs  from 
the  time  required  for  wave  propagation  through  the  substructure. 

If  full  cancellation  of  incoming  sinusoidal  disturbances  is 
to  be  achieved  for  the  portion  of  the  substructure  beyond  the 
controller,  the  controller  feedforeward  gain  amplitude  must  match 
that  of  the  substructure  exactly.  If  attenuation  is  tc  be 
achieved,  the  controller  gain  amplitude  must  be  less  than  twice 
that  of  the  substructure. 

For  a  nonzero  discrepancy  between  the  controller  delay  and 
the  delay  for  wave  propagation  in  the  structure,  the  ratio  cf  the 
system  (sinusoidal)  output  amplitude  to  the  input  sinusoidal 
amplitude  is  dependent  on  the  input  waveform's  radian  frequency. 

To  fully  cancel  sinsuoidal  disturbances  in  the  substructure 
beyond  the  controller,  assuming  matched  controller  and 
substructure  gains,  the  duration  of  the  controller-substructure 
delay  discrepancy  must  be  equal  to  an  integer  multiple  of  naif 
the  input  waveform  period.  To  attenuate  incoming  disturbances, 
assuming  matched  controller  and  substructure  gains,  the  delay 


discrepancy  must  not  differ  from  an  integer  multiple  of  the 


input  waveform  period  by  more  than  one  sixth  of  the  period 


duration.  This  range  approaches  zero  as  the  magnitude  of  K 


approacnes  2,  and  becomes  larger  as  the  magnitude  of  K  approaches 


Further  work  should  be  done  to  determine  analytically  the 


dependence  on  frequency  and  on  controller  error  of  the  ratio  of 


output  to  input  wave  energies.  The  effects  of  delay  error  in  the 


stabilizing  feedback  loop  should  be  analyzed. 


The  analytical  results  contained  here  should  also  be 


verified  experimentally. 
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Fig.  4  Block  diagram  with  superposition  at  cross  sections. 
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APPENDIX  A 


ANALYSES  OF  FEEDBACK  CONTROL  OF  WAVE  PROPAGATION  IN  STRUCTURES 

The  feedback  control  of  wave  propagation  in  lattice 
structures  is  considered  in  this  appendix.  Limitations  of  such  a 
control  scheme  are  illustrated  and  discussed. 

The  geometry  of  a  one-dimensional  lattice  substructure, 
which  is  modelled  as  a  thin  rod,  is  shown  in  Fig.  A-1 .  As  shown, 
the  coordinate  x  refers  to  the  undisturbed  location  of  a  cross 
section  of  the  substructure.  Two  sections  of  the  substructure, 

separated  by  distance  1  ,  are  labelled  as  x  and  x  .  A 

12  12 

section  to  the  right  of  x  is  labelled  as  x  .  The  section  x 

2  3  3 

represents  the  location  in  the  substructure  to  be  protected  from 
incoming  stress  disturbances  by  the  feedback  Bystem. 

In  the  feedback  control  scheme,  a  disturbance  is  assumed  to 
originate  from  the  left  of  x  and  propagate  through  the 
substructure.  A  sensor  measures  the  stress  (or  displacement  or 
velocity)  at  cross  section  xz ,  and  the  feedback  system  provides  a 
force  in  the  longitudinal  direction  at  cross  section  x^  through 
an  actuator.  The  force  provided  by  the  actuator  generates  a 
rightwardly-propagating  stress  wave  at  x  ^  whose  instantaneous 

amplitude  at  x  is  the  exact  negative  of  the  stress  sensed  at  x  . 

i  2 

The  force  provided  by  the  actuator  also  generates  a  corresponding 
lef twardly-propagating  stress  wave  at  x  which  is  the  negative  of 
the  rightwardly-propagating  stress  wave.  (In  this  study, 
lef twardly-propagating  waves  are  not  reflected  at  the  left 
boundary  of  the  substructure.)  The  total  response  of  the 
substructure  is  the  superposition  of  the  incoming  stress  wave 
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disturbance  and  the  stress  waves  generated  by  the  actuator. 

Figs.  A-2  and  A-3  show  the  responses  of  the  substructure 
under  feedback  control  to  an  input  single  pulse  disturbance  and 
to  an  input  step  disturbance,  respectively.  The  sensor  S  is 

located  at  x  and  the  actuator  A  is  located  at  x  .  The  feedback 

2  1 

system  is  indicated  by  a  block  diagram  containing  -1  as  the  only 
entry.  The  transit  time  for  wave  propagation  from  x;  to  x  is 

x  .  Figs.  A-2  and  A-3  show  the  stress  distributions  a  (x,t) 

1  2 

along  the  length  of  the  substructure  at  equal  time  intervals  for 
increasing  time.  The  time  interval  shown  corresponds  to  half  the 

transit  time  t  A  rightwardly-propagating  stress  wave 

12 

disturbance  originates  to  the  left  of  x:  •  The  stress 
distribution  along  the  substructure  is  observed  at  tine  t  , 
after  the  disturbance  has  passed  location  x  but  before  it  has 

l 

reached  location  x  .  Then,  the  stress  distributions  along  the 

2 

substructure  att,t,...,t,t  are  shown. 

1  2  6  7 

In  Figs.  A-2  and  A-3 »  the  solid  outlines  represent  the 
incoming  stress  disturbance  and  the  dashed  outlines  represent  the 
stress  waves  generated  by  the  actuator.  Also,  the  shaded 
portions  represent  the  net,  or  physically  realized  magnitude  of 
the  stress  waves  in  the  substructure  due  to  superposition  of  the 
incoming  disturbance  and  the  actuator-generated  stress  waves. 

Figs.  A-2a  through  A-2h  show  the  responses  of  the 
substructure  to  a  single  pulse  disturbance.  The  duration  of  the 

disturbance  is  assumed  to  be  less  than  the  transit  time  t  .  Fig. 

12 

A-2b  shows  the  incoming  disturbance  just  before  it  arrives  at  the 
sensor  S.  Fig.  A-2c  shows  the  incoming  disturbance  as  it  is  just 
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passing  S  and  the  actuator  is  active  in  producing  both  a 
rightwardly-propagating  wave  and  a  lef twardly-propagating  wave. 
Fig.  A-2d  shows  the  incoming  disturbance  after  it  has  passed  S 
and  the  actuator  is  not  active.  However,  a  negative  pulse 
propagating  to  the  right  has  been  generated  by  the  actuator. 
Fig.  A-2e  shows  the  actuator  generating  stress  waves  due  to  the 
negative  actuator-generated  rightwardly-propagating  wave.  This 
process  continues  indefinitely  and  an  infinite  sequence  of 
identical  pulses,  with  alternating  sign  and  separated  by  T  ,  is 

sent  to  the  right  of  location  x  .  Thus,  if  a  point  to  the  right 

2 

of  x  (that  is,  point  x  )  is  to  be  protected  from  stress 

2  3 

disturbances  by  the  feedback  system,  the  feedback  system  actually 
worsens  the  situation  by  generating  more  pulses  towards  the  point 
of  interest. 

Figs.  A-3a  through  A-3h  show  the  responses  of  the 
substructure  to  a  step  disturbance.  Fig.  A-3a  shows  the  incoming 
disturbance  just  before  it  arrives  at  the  sensor  S.  Fig.  A-3b 
shows  the  leading  edge  of  the  incoming  disturbance  has  just 
passed  S  and  the  actuator  is  active  in  producing  both  a 
rightwardly-propagating  wave  and  a  lef twardly-propagating  wave. 
Tne  actuator-generated  rightwardly-propagating  wave  cancels  the 
portion  of  the  incoming  disturbance  that  overlaps  with  it.  Fig. 
A-3a  shows  the  leading  edge  of  the  actuator-generated 
rightwardly-propagating  wave  has  passed  the  sensor  S  and  has 
cancelled  the  incoming  disturbance  at  S.  Eecause  there  is  no 
longer  any  stress  at  S,  the  actuator  is  not  active.  Thus,  the 
incoming  disturbance  arriving  at  the  actuator  A  passes  unaltered 
as  shown  in  Figs.  A-3d  and  A-3e. 


Fig.  A-3f  shows  the  leading 


eage  of  the  unaltered  incoming  disturbance  has  just  passed  the 
sensor  S  and  the  actuator  is  active  again.  This  process 

continues  indefinitely  and  an  infinite  sequence  of  identical 

pulses,  each  with  duration  equal  to  the  transit  time  t  and  each 

22 

also  separated  by  t  ,  is  sent  to  the  right  of  location  x 

12  2 

Thus,  if  a  point  to  the  right  of  x  (that  is,  point  x  )  is  to  be 

2  3 

protected  from  stress  disturbances  by  the  feedback  system,  the 
feedback  system  does  not  remove  all  of  the  step  stress 

disturbance,  but  does  prevent  segments  of  the  input  step  from 
reaching  the  point  of  interest.  The  penalty  is  that  segments  of 
lef twardly-propagating  waves  to  the  left  of  x  ^  have  double  the 
amplitude  of  the  incoming  disturbance. 

To  summarize,  the  feedback  control  of  wave  propagation  in 
lattice  substructures  as  described  here  has  significant 
limitations.  For  short  input  pulse  disturbances  as  shown  in  Fig. 
A-2,  an  infinite  sequence  of  pulses  of  alternating  sign  is 

generated  to  the  right  of  location  x  .  For  long  input  pulse 

2 

disturbances  as  shown  in  Fig.  A-3>  only  a  portion  of  the 
disturbance  is  prevented  from  reaching  a  location  to  the  right  of 
x  (that  is,  point  x  ).  However,  as  shown  in  Figs.  A-2  and  A-3, 

2  3 

this  control  system  does  not  appear  to  cause  system  instability 
because  the  stress  amplitudes  are  always  limited.  Nevertheless, 
feedback  control  may  still  be  useful  as  a  part  of  a  total  control 
scheme  for  the  control  of  wave  propagation  in  lattice  structures. 


APPENDIX  B 


IMPULSE  RESPONSE  AND  TRANSFER  FUNCTION  FOR 
WAVE  PROPAGATION  FEEDFORWARD  CONTROL 

The  impulse  response  and  the  transfer  function  for  the 
feedforward  control  of  wave  propagation  in  structures  are 
considered  in  this  Appendix. 

In  general,  a  linear  system  can  be  characterized  by  its 
impulse  response  h(t)  or  its  transfer  function  H(u>)  as  shown  in 
Fig.  B-1  [4].  (The  transfer  function  is  also  called  the  frequency 
response.)  For  an  input  x(t),  where  t  denotes  time,  the  output 
y(t)  can  be  obtained  via  h(t)  as  [4] 

y(t)  =  f  x(T)h(t-x)dT  (B.l) 

—  OO 

where  T  is  a  variable  of  integration.  For  an  input  X(w),  where 
w  denotes  radian  frequency,  the  output  Y(w)  can  be  obtained  via 
H (oj  )  as  [4] 

Y(w)  =  H(oj)  X(w)  (B. 2) 


Furthermore,  the  impulse  response  and  the  transfer  function  are 
related  as  [4] 


H(w) 


f  h(t> 


(B.3) 


where  i  denotes  /- 1  . 


Fig.  B-2  shows  a  schematic  of  the  feedforward  control  of 


wave  propagation  in  structures.  The  substructure  is 


characterized  by  a  time  delay  of  t 


i  2 


The  feedforward 


controller  is  characterized  by  a  time  delay  of  t  '  and  an 

12 

amplification  factor  K.  (Waves  may  be  generated  at  the  actuator 
and  propagate  back  along  the  structure  toward  the  sensor. 
However,  it  is  assumed  in  this  Appendix  that  these  particular 
waves  are  ignored  by  the  sensor.  Also,  it  is  noted  that 
dispersion  and  attenuation  are  not  considered  in  the  model  shown 
in  Fig.  B-2.)  The  output  of  the  feedforward  control  system  is 
formed  by  subtracting  the  output  of  the  feedforward  controller 
from  the  output  of  the  substructure  as  indicated  in  Fig.  B-2. 

Considering  an  input  impulse  6(t)  into  the  feedforward 
control  system  shown  in  Fig.  B-2,  the  impulse  response  hf(t)  of 
the  feedforward  control  system  can  be  written  as 


h,(t)  =  6 ( t—T  )  -  K<5 ( t- T '  ) 
t  12  12 


(B.4) 


where  5  is  the  Dirac  delta  function  [4]. 

Eqn.  (B.4)  shows  that  the  output  of  the  feedforward  control 
system  is  formed  by  subtracting  the  output  of  the  feedforward 

controller,  which  is  the  input  impulse  time-delayed  by  x  ’  and 

1 2 

scaled  by  K,  from  the  output  of  the  substructure  only,  which  is 

the  input  impulse  time-delayed  by  x  and  scaled  by  unity. 

1 2 

Then,  the  transfer  function  H^(u))  of  the  feedforward  control 
system  shown  in  Fig.  B-2  can  be  obtained  from  h^.(t)  via 
eqn.  (E.3)  as 
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Hf(oi)  = 


(B.5) 


Substituting  eqn.  (B.4)  into  eqn .  (£.5)  gives 


Hf(u>) 


co 

f  [6(t-T12)-K6(t-TipJ  e-iWtdt 


(B 


From  the  properties  of  Dirac  delta  function,  eqn.  (B.6)  can 
integrated  to  give  [4] 


Hf  (w) 


■  1C0T 


1  2 


(B 


Eqn.  (B.7)  gives  the  transfer  function  Hf( w  )  of 
feedforward  control  system  shown  in  Fig.  B-2.  As  shown 
eqn.  (B.7),  the  transfer  funtion  is  a  complex  function, 
magnitude  of  the  transfer  function  will  be  discussed. 

The  square  of  the  magnitude  of  Hf(w)  can  be  obtained 
[4] 

|  Hf  (tu)  |  2  =  Hf(w)Hf(-w)  (B 


where  |  H^(u)|  denotes  the  square  of  the  magnitude  of  Hf(w) 
Hf(- w)  denotes  the  complex  conjugate  of  Hf(w). 

The  complex  conjugate  of  Hf(w)  can  be  obtained  from 
(B.7)  as 


Hf  (— cj) 


1WT,  „ 
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Substituting  eqns.  (B.7)  and  (B.9)  into  eqn.  (B.8)  gives 


4 


Eqn.  (£.10)  can  be  expanded  as 


| Hf  (co)  |  2  =  e_iwT  j 2  +  Euit  j 2  _  Ke-iu>xi2  +  iuV2 


-  Ke~iU)Ti2  +  iwTi2  +  K2e"iWTi2  +  ia)Ti2 


(B. 11) 


Eqn.  (£.11)  can  be  simplified  as 


| Hf  (to)  |  2  =  l-Keia)(Ti2-Ti2)  -  Ke~1U)(Ti2-Ti2)  +  K2 


(B . 12) 


The  second  and  third  terms  on  the  right-hand  side  of  eqn.  (£.12) 
are  complex  conjugates  of  each  other.  When  added,  the  imaginary 
parts  of  these  two  terms  cancel.  Thus,  eqn.  (£.12)  becomes 


| Hf  (u))  j  2  =  1  -2Kcos  T12)]  +  K2 


(B.13) 


Defining 


t  '  -  t  =  e 
12  12 


(B.14) 


where  £  denotes  the  difference  in  time  delays  between  the 
feedforward  controller  and  the  substructure,  eqn.  (£.13)  can  be 
rewritten  as 


| Hf  (cu)  |  =l-2Kcos(we)  +  K2 


(B. 15) 


Eqn.  ( B . 1 5 )  gives  the  square  of  the  magnitude  of  the 
transfer  function  for  feedforward  control  of  wave  propagation  in 


K.  -J 

•\V>> 


L  *  .N 

•  .  •  o  •  * 


•  A 

•  .  ’»  .  ** 

^ySv- 

**.**•■«»  ** 
•>  -  -V 

■**_ 

•  • 

t 


%  ~4 


a  .• 
w'aV’V 


AVY; 

the 

A  • -  *% 

in 

«  .  «<  •  «r\. 

*  * .  .  *  w 

*‘avy 

An'.'*, 


.  a  a 

•, 

V/fcN\* 


the  simple  structures  under  consideration.  Eqn.  (B.15)  is 
plotted  in  Fig.  B-3  versus  the  parameter  (ue )  . 

To  summarize,  the  eqn.  (B.15)  and  Fig.  B-3  account  for  the 
general  case  where  the  time  delay  in  the  feedforward  controller 
exceeds  the  time  delays  in  the  substructure  by  e  and  also  where 
the  feedforward  amplification  is  K,  as  shown  in  Fig.  B-2. 
Fig.  B-3  shows  that  the  magnitude  of  the  transfer  function  is 
periodic  in  we  .  Thus,  when  we  is  increased  (or  decreased)  by 
2v  ,  the  transfer  function  magnitude  repeats  itself.  Also,  for 
a  given  u  ,  a  positive  e  and  a  negative  £  of  the  same  magnitude 
produce  the  same  transfer  function  magnitude.  Furthermore,  Fig. 
E-3  shows  that  for  any  value  of  £  the  magnitude  of  the  transfer 
function  has  the  least  possible  value  when  K  equals  unity. 
Finally,  because  the  transfer  function  magnitude  as  shown  in  Fig. 
b-3  is  frequency  dependent,  the  feedforward  control  system 
response  due  to  certain  frequencies  (for  constant  £  )  is  enhanced 
wnile  tnat  from  other  frequencies  is  diminished. 


Eqn.  ( B.1 0 )  can  be  expanded  as 


| Hf  (oj)  |  2  =  e-iUTi2+iwTi2  -  + 


-  Ke“iU)Ti2  +  ia>Ti2  +  K2e_iU)Ti2  +  iU)Ti2 


(B.ll) 


Eqn.  ( B. 1 1 )  can  be  simplified  as 

| Hf  (oo)  |  2  =  l-Keia)(Ti2_Ti2)  -  Ke~iW^Ti2~Ti2^  +  K2  (B.12) 

The  second  and  third  terms  on  the  right-hand  side  of  eqn.  (B.12) 
are  complex  conjugates  of  each  other.  When  added,  the  imaginary 
parts  of  these  two  terms  cancel.  Thus,  eqn.  (B.12)  becomes 


|Hf(u) 


-2Kcos 


+  K 


(B.13) 


Defining 

x  '  -  t  =  e  (B.14) 

12  12 

where  £  denotes  the  difference  in  time  delays  between  the 
feedforward  controller  and  the  substructure,  eqn.  (B.1J)  can  be 
rewritten  as 


| (to)  ]  =  1  -2Kcos(coe)  +  K2  (B.15) 

Eqn.  (B.15)  gives  the  square  of  the  magnitude  of  the 
transfer  function  for  feedforward  control  of  wave  propagation  in 


the  simple  structures  under  consideration.  Eqn.  (£.15)  Is 
plotted  in  Fig.  B-3  versus  the  parameter  (uae  )  . 

To  summarize,  the  eqn.  ( B . 1 5 )  and  Fig.  B-3  account  for  the 
general  case  where  the  time  delay  in  the  feedforward  controller 
exceeds  the  time  delays  in  the  substructure  by  e  and  also  where 
the  feedforward  amplification  is  K,  as  shown  in  Fig.  B-2 . 
Fig.  B-3  shows  that  the  magnitude  of  the  transfer  function  is 
periodic  in  we  .  Thus,  when  we  is  increased  (or  decreased)  by 
2tt  ,  the  transfer  function  magnitude  repeats  itself.  Also,  for 
a  given  w  ,  a  positive  e  and  a  negative  £  of  the  same  magnitude 
produce  the  same  transfer  function  magnitude.  Furthermore,  Fig. 
B-3  shows  that  for  any  value  of  £  the  magnitude  of  the  transfer 
function  has  the  least  possible  value  when  K  equals  unity. 
Finally,  because  the  transfer  function  magnitude  as  shown  in  Fig. 
B-3  is  frequency  dependent,  the  feedforward  control  system 
response  due  to  certain  frequencies  (for  constant  £  )  is  enhanced 
while  that  from  other  frequencies  is  diminished. 


